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Role  of  Geometry  on  the  Color  of 
Flux  Noise  in  dc  SQUIDs 

F.  C.  Wellstood,  C.  Urbina,  and  John  Clarke 


Abstract — We  examine  the  behavior  of  low  frequency  flux  noise 
measured  in  ten  Nb  —  XbO„  —  Pbln  dc  SQUIDs  in  the  temperature 
range  1  K  to  20  mK.  As  previously  reported,  the  flux  noise  power 
spectrum  S$(f)  typically  scales  with  frequency  as  A/f".  Remark¬ 
ably,  the  excess  noise  power  A  at  1  Hz  increases  as  the  tempera¬ 
ture  is  decreased  below  about  0.7  K,  saturating  to  a  value  around 
(10  ,($0)2/Hz  below  about  0.2  K.  Here  we  report  on  the  depen¬ 
dence  of  both  the  magnitude  A  and  the  slope  a  on  the  size  and 
linewidth  of  the  SQUID  loops.  In  particular,  at  the  lowest  temper¬ 
atures  we  find  that  a  declines  to  values  as  low  as  0.5  in  the  smallest 
devices,  while  a  is  in  the  range  0.9-1  for  the  largest  devices. 

Index  Terms — Flux  noise,  SQUIDs,  two-level  systems,  1/f  noise. 


I.  Introduction 

RECENT  measurements  on  flux  qubits  reveal  the  presence 
of  excess  low-frequency  flux  noise  that  causes  dephasing 
when  the  device  is  biased  away  from  the  “sweet-spot”  [1].  This 
flux  noise  appears  to  be  of  the  same  type  as  that  first  reported 
more  than  20  years  ago  in  a  number  of  thin-film  dc  SQUIDs  that 
were  cooled  below  1  K  [2],  [3].  The  spectral  density  S$(f)  of 
this  noise  was  found  to  be  of  the  form  A/f“,  where  f  is  the  fre¬ 
quency  and  ?  typically  lies  between  0.5  and  1.0.  Furthermore, 
A  increased  rather  slowly  with  the  dimensions  of  the  SQUIDs. 
Only  recently — with  the  impetus  of  research  into  quantum  com¬ 
puting — has  significant  new  theoretical  and  experimental  work 
been  undertaken  to  identify  the  underlying  cause  of  the  noise. 
A  consensus  has  been  developed  that  the  noise  is  due  to  ran¬ 
domly  flipping  electronic  spins  distributed  in  a  thin  layer  on  the 
metallic  surfaces  of  the  devices  [4] — [13].  Theoretical  investiga¬ 
tions  [4],  [9],  [12]  showed  that  this  picture  is  consistent  with 
the  slow  scaling  of  the  noise  power  with  SQUID  dimensions. 
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In  particular  [9],  [12],  it  was  shown  that  S$(f)  scales  as  d/w  in 
the  limit  w/d  -C  1 ,  where  d  is  the  outer  dimension  of  the  SQUID 
and  w  is  the  linewidth  of  the  SQUID  loop.  This  scaling  was  con¬ 
firmed  experimentally  by  Lanting  et  al.  [13]. 

Despite  this  progress,  many  unusual  features  of  the  noise  re¬ 
main  unexplained  or  poorly  understood.  Here,  we  examine  addi¬ 
tional  results  from  our  early  study  [2]  not  previously  published 
that  highlight  some  of  the  most  puzzling  features  of  flux  noise,  in 
particular  the  unusual  dependence  of  the  exponent  a  on  the  tem¬ 
perature  and  geometry  of  the  device.  In  Section  II  we  describe 
the  SQUIDs  that  we  studied,  and  in  Section  III  briefly  describe 
our  measurement  technique.  We  present  our  experimental  re¬ 
sults  in  Section  IV,  and  in  Section  V  discuss  them  in  terms  of 
a  loss  tangent.  We  describe  the  temperature  and  geometry  de¬ 
pendence  of  a  in  Section  VI,  and  present  our  conclusions  in 
Section  VII. 

Drung  et  al.  [14]  reported  similar  variations  in  a  at  the  same 
Applied  Superconductivity  Conference. 

II.  SQUID  Geometry  and  Parameters 

To  understand  how  the  flux  noise  varies  with  the  SQUID 
layout,  we  fabricated  dc  SQUIDs  with  a  variety  of  sizes  and 
geometries.  Figs.  1  and  2  show  the  layout  of  the  larger  and 
smaller  SQUIDs,  respectively.  The  SQUIDs  were  fabricated  on 
oxidized  Si  wafers  (except  for  devices  K1  and  J1  which  were 
on  a  sapphire  wafer)  using  photolithographic  patterning.  First, 
a  35-nm  thick  AuCu  resistive  shunt  layer  was  evaporated  and 
patterned  with  liftoff.  This  was  followed  by  a  200-nm  thick  sput¬ 
tered  Nb  layer  which  was  patterned  in  an  SF6  +  O2  plasma  etch. 
Two  SiO  insulation  layers  (500-nm  total  thickness),  were  evap¬ 
orated,  using  a  1-2  nm  underlayer  of  Cr  or  Ti  as  an  adhesion 
promoter,  and  2x2  (pm)2  windows  for  the  tunnel  junctions 
were  opened  using  liftoff.  The  counter-electrode  layer  was  pat¬ 
terned  in  resist  and  the  wafer  scribed  into  chips. 

The  tunnel  junctions  were  formed  by  placing  a  chip  in  a  dif¬ 
fusion-pumped  vacuum  chamber  and  ion-milling  the  exposed 
Nb  electrode  for  40  s  to  remove  surface  oxide.  The  chamber 
was  vented,  and  the  chip  was  removed  and  remounted  in  an 
rf  Ar  —  O2  plasma  oxidation  chamber  to  produce  NbOx  on 
the  cleaned  Nb  electrode.  The  chamber  was  re-evacuated  and 
a  200-nm  layer  of  Pb+5%In  evaporated  to  form  the  junction 
counter-electrodes,  which  were  patterned  with  lift-off  in  ace¬ 
tone. 

Table  I  gives  the  parameters  for  ten  devices.  The  loop  induc¬ 
tance  F  was  estimated  from  the  measured  modulation  depth  of 
the  current  I  versus  flux  <!>;  this  introduced  some  spread  in  the 
values,  even  for  devices  with  identical  geometry.  For  example, 
although  Dl,  D2,  Ml,  and  M2  have  the  same  loop  geometry. 
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Fig.  1.  (a)  Side-view  of  layer  structure  used  in  the  SQUIDs.  A  thin  AuCu  layer 
is  deposited  on  a  substrate  and  patterned  to  form  shunts  (not  shown).  This  is 
followed  by  a  Nb  layer  and  two  insulating  layers  of  SiO  that  define  two  2  /mi  x 
2  /mi  windows.  A  NbOx  tunnel  barrier  is  grown  on  the  Nb  and  the  top  layer 
of  Pb  +  5%  In  is  added  to  complete  the  Nb  —  NbOx  —  Pbln  tunnel  junctions. 
(b)-(d)  Configuration  of  larger  SQUIDs:  Types  D,  K,  and  P.  Type  M  SQUIDs 
have  the  same  layout  as  type  D,  except  for  the  addition  of  large  area  AuCu 
cooling  fins  to  the  shunts  [15]. 
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Fig.  2.  Configuration  of  smaller  SQUIDs:  Types  C,  F,  J  and  O.  Junction  and 
layer  descriptions  same  as  in  Fig.  1 . 


TABLE  I 

Device  Parameters 


# 

L 

2Io 

d 

w 

i 

Ap 

-Ja 

a 

(pH) 

(pa) 

(pm) 

(pm) 

(pm) 

(pm)3 

•JJk 

D1 

440 

6.3 

1000 

300-550 

2200 

2x1  O' 

14.5 

0.96 

D2 

500 

4.5 

1000 

300-550 

2200 

2x1 0s 

14.5 

0.95 

Ml 

515 

5.6 

1000 

300-550 

2200 

2x1 05 

18 

0.94 

M2 

470 

6.2 

1000 

300-550 

2200 

2xl05 

14 

0.92 

PI 

1500 

1.77 

144-740 

10 

2560 

5xl04 

14.5 

0.74 

K.1 

530 

11.4 

75-1000 

20 

1700 

7xl04 

11 

0.72 

.11 

200 

21.6 

200 

50 

560 

2xl04 

5 

0.71 

FI 

700 

5.7 

100 

2 

390 

IxlO4 

8 

0.67 

01 

45 

1 .36 

90-105 

36 

120 

2xl03 

7 

0.60 

C2 

210 

5.3 

200 

50 

400 

2xl04 

7.5 

0.50 

SQUID  parameters  listed  in  descending  order  of  the  slope  a .  First  column 
gives  the  device  name.  L  is  the  loop  inductance,  2I0  is  the  maximum 
critical  current  (for  nominally  zero  flux  bias),  d  is  the  outside  dimension  of 
the  SQUID  loop  (two  dimensions  are  given  for  non-square  loops),  w  is  the 
linewidth  of  the  SQUID  loop,  f  is  the  approximate  length  around  the  inner 
hole  of  the  SQUID  loop,  and  Ap  is  the  approximate  pick-up  area  of  the 
loopVd  is  the  square  root  of  the  excess  flux  noise  power  spectral  density  at 
1  Hz  and  a  is  the  noise  slope,  both  taken  at  the  lowest  temperature  measured 
for  the  device. 


their  estimated  loop  inductances  range  from  440  pH  to  515  pH 
with  an  average  of  480  pH. 

III.  Experimental  Arrangements 

The  low-frequency  magnetic  flux  noise  signal  from  a  dc 
SQUID  is  small.  To  measure  this  noise  we  connected  the 
device  under  test  (SQUID  1)  in  series  with  a  small  resistor 
Rx  ~  0.07  D,  an  rf-choke  (Lrf  ~  12  //Hj  and  the  input  coil 
Li  of  a  second  dc  SQUID  acting  as  a  noise  detector  (Fig.  3). 
SQUID  2  was  a  Nb  square  washer  design  with  1  mm  outside 
dimension,  a  180  /im  inner  hole  dimension  and  an  integrated 
input  coil  with  20  turns  (except  for  that  used  for  measuring 
device  M2  which  had  50  turns).  This  SQUID  was  operated 
in  a  flux-locked  loop  (FLL)  [16],  a  current  source  Ri  was 
connected  across  Rx,  and  a  second  current  source  I,\,\  was  used 
to  apply  flux  to  SQUID  1 . 

For  a  given  applied  flux  <T>  | ,  SQUID  1  was  voltage  biased  by 
applying  current  Ri  greater  than  its  critical  current.  Changes  in 
flux  in  SQUID  1  change  the  current  R  passing  through  it  and 
Li,  causing  the  flux  <E>o  in  SQUID  2  to  change  and  ultimately 
producing  a  change  in  the  output  voltage  Vout  from  the  FLL. 
The  arrangement  allows  us  to  measure  the  noise  in  SQUID  1  as 
well  as  its  current-voltage  (R  —  VR  and  current-flux  (R  —  <LR 
characteristics. 

To  obtain  the  power  spectral  density  of  the  system  noise  at 
specific  R  —  'R  bias  points  for  SQUID  1,  we  connected  the 


Fig.  3.  Flux  noise  in  SQUID  1  is  measured  by  SQUID  2  operated  in  a  flux- 
locked  loop.  SQUID  1  is  voltage-biased  by  Rx  =  0.07  Q  and  current  source 
I^i  so  that  changes  in  the  flux  <[>!  in  SQUID  1  lead  to  changes  in  the  current  Ix 
through  SQUID  1.  Changes  in  Vout  from  the  FLL  are  proportional  to  changes 
in  flux  in  SQUID  2,  which  in  turn  are  proportional  to  changes  in  Ii. 


voltage  output  of  the  flux-locked  loop  to  a  spectrum  analyser. 
From  the  slope  of  the  R  —  <L|  curve  at  the  bias  point,  this  noise 
power  spectrum  is  readily  transformed  into  an  effective  flux 
noise  in  SQUID  1  (Fig.  4). 

The  measured  and  measurement  SQUIDs  were  enclosed  in 
separate  compartments  in  a  Nb  tube  mounted  in  a  sealed  Cu 
cell.  The  cell  was  bolted  to  the  mixing  chamber  of  a  dilution 
refrigerator  and  filled  with  superfluid  helium  for  good  thermal 
contact.  To  reduce  the  ambient  magnetic  field,  external  magnetic 
noise  and  rf-interference,  the  dewar  was  surrounded  with  two 
mu-metal  cylinders  and  the  entire  experiment  was  installed  in  a 
Cu-mesh  rf-shielded  room.  We  made  noise  measurements  from 
4.2  to  0.02  K,  with  SQUID  1  biased  to  a  few  //,V  (power  ~ 
0.1  nW). 

IV.  Flux  Noise  Spectral  Density  Versus  Temperature 

Fig.  4  shows  two  typical  flux  noise  power  spectra.  Ignoring 
noise  peaks  at  60  and  180  Hz,  the  spectra  obey 

S*(/)  =  A  f~a  +  C, 


(1) 
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f(Hz) 

Fig.  4.  Power  spectral  density  of  magnetic  flux  noise  S<$  versus  frequency  f 
for  SQUID  FI  (open  triangles)  at  T  =  0.14  K  and  M2  (closed  circles)  at  T  = 
0.021  K.  Dashed  curve  is  fit  to  white  noise  plus  1  /f“  noise  with  a  =  0.67  for 
FI  and  solid  curve  is  similar  fit  for  M2  with  a  =  0.92. 


Fig.  5.  Magnitude  of  magnetic  flux  noise  at  1  Hz  versus  temperature  T.  The 
symbols  are  for  devices  D2  -BrMl  -^-Dl  -*tM2  -tvPl  —’fc-Ji- 

— ©--FI-— A,-Oi  -*randC2  -*rr 

where  A  is  a  constant,  C  is  frequency-independent  (white)  noise 
from  the  device  under  test  and  measuring  system,  and  a  is  the 
slope  of  the  excess  noise  on  a  log-log  plot.  In  analogy  with  white 
noise,  noise  with  a  near  1  is  called  “pink  noise”  because  there 
is  more  noise  at  lower  (“red”)  frequencies.  As  Fig.  4  shows,  the 
noise  in  SQUID  M2  is  more  “red”  than  that  in  FI. 

Fig.  5  shows  the  square  root  of  the  flux  noise  power  spectrum 
at  1  Hz  plotted  versus  temperature  T  for  the  10  devices  listed  in 
Table  I.  For  T  >  1  K,  the  noise  may  contain  significant  con¬ 
tributions  from  critical  current  noise  [17],  but  this  falls  rapidly 
as  the  temperature  decreases.  Below  1  K  the  excess  noise  was 
dominated  by  a  flux-like  term,  that  is,  the  output  noise  scaled 
with  the  slope  of  the  I  |  —  <f>i  curve,  as  expected  for  a  fluctuating 
flux  in  the  SQUID. 

As  noted  previously,  there  are  remarkable  features  in  this  be¬ 
havior.  (i)  Broadly  speaking,  devices  with  similar  geometry  tend 
to  produce  similar  noise,  (ii)  Ignoring  the  behavior  above  1  K 
which  contains  significant  critical  current  noise  [17]  and  is  ma¬ 
terial  dependent  [2],  the  noise  tends  to  increase  as  T  decreases 
from  1  K.  (iii)  The  largest  devices  (Dl,  D2,  Ml,  M2,  bolder 
solid  curves  in  Fig.  5)  tend  to  have  larger  noise  with  abroad  peak 
around  0.5  to  0.7  K.  Smaller  devices  (for  example,  C2  and  01, 
thin  solid  curves  in  Fig.  5)  have  smaller  noise,  which  tends  not 
to  change  much  with  temperature  below  1  K  or  increase  steadily 
as  T  decreases,  (iv)  At  the  lowest  temperatures  the  noise  for  all 
the  devices  falls  in  the  range  5  to  18  //'ho/H'/1,2  at  1  Hz,  de¬ 
spite  the  diverse  range  of  device  parameters  (Table  I):  the  in¬ 
ductances  L  range  from  45-1500  pH,  the  critical  currents  range 
from  1-22  //  A,  the  pickup  areas  Ap  of  the  SQUID  loops  vary 


from  about  2  X  103  (pm)2  to  2  X  105  (pm)2,  and  the  linewidths 
of  the  films  used  to  make  the  SQUID  loops  vary  from  2  pm  to 
300  pm  or  greater.  In  addition,  devices  (J1  and  Kl)  made  on 
a  sapphire  substrate  also  exhibit  similar  noise,  and  the  fraction 
of  the  loop  made  from  Nb  or  Pbln  does  not  seem  to  affect  the 
noise  level.  These  general  observations  are  consistent  with  mea¬ 
surements  on  other  devices,  embracing  a  wider  range  of  param¬ 
eters,  including  measurements  by  Yoshihara  et  al.  [1]  on  much 
smaller  aluminum  flux  qubits  [Ap  ~  2  (pm)2]  with  somewhat 
smaller  estimated  levels  of  noise  (1-2  //'ho/Hz1  ,/2  at  1  Hz).  As 
mentioned  earlier,  the  relatively  slow  dependence  of  the  noise 
amplitude  on  the  SQUID  dimensions  can  be  explained  in  terms 
of  fluctuating  electron  spins  distributed  uniformly  on  the  surface 
[4],  [9],  [12], 

V.  Magnetic  Loss  Tangent 

The  existence  of  spin  fluctuations  must  be  accompanied  by 
dissipation.  For  a  magnetic  material,  we  can  characterize  the 
dissipation  in  terms  of  the  magnetic  loss  tangent,  tan(pm),  the 
ratio  of  the  imaginary  part  of  the  magnetic  susceptibility  to  the 
real  part.  Applying  the  fluctuation-dissipation  theorem  to  a  ma¬ 
terial  with  a  magnetic  loss  tangent  yields: 

tan(<5m)  ss  S$(f)/2kBTLrj.  (2) 

The  filling  factor  //  accounts  for  the  distribution  of  material 
around  the  SQUID  loop:  //  =  I  for  a  uniform  3D  distribution. 

Eq.  (2)  implies  that  tan(pm  )  is  independent  of  frequency  if 
the  noise  scales  as  1/f  and  that  loops  with  the  same  filling  factor 
will  exhibit  a  noise  power  spectrum  that  scales  with  the  loop 
inductance.  We  can  view  tan(pm  )  as  a  figure  of  merit  that  ac¬ 
counts  for  expected  effects  of  the  loop  and  source  geometry. 
Here,  as  the  distribution  of  the  sources  is  unknown,  we  set  r\  —  1 
to  find  an  effective  tan(pm). 

Fig.  6  is  aplot  of  the  effective  tan(pm)  at  1  Hz  versus  temper¬ 
ature  T  for  the  noise  data  shown  in  Fig.  5.  Overall,  the  magnitude 
is  rather  low,  but  increases  significantly  as  the  temperature  de¬ 
creases.  Fig.  6  shows  that  at  any  temperature  the  loss  tangent  of 
the  devices  falls  within  a  band  that  spans  an  order  of  magnitude. 
While  both  Figs.  5  and  6  show  a  spread  of  an  order  of  magni¬ 
tude,  note  that  Fig.  5  is  a  plot  of  the  square  root  of  the  power 
spectral  density,  while  tan (/im)  is  proportional  to  the  power. 
Although  further  analysis  of  the  filling  factor  is  needed,  this  be¬ 
havior  implies  that  the  noise  power  scales  approximately  with 
loop  inductance  L.  Thus,  the  noise  energy  S$(f) /2  L  is  roughly 
constant  at  a  given  temperature. 

VI.  Slope  of  Flux  Noise  Versus  Temperature 

Fig.  7  shows  aplot  of  the  slope  a  of  the  1/f"  flux  noise  versus 
temperature  T  for  the  10  devices.  The  values  of  a  at  the  lowest 
temperatures  are  also  listed  in  Table  I.  The  four  largest  devices 
(Dl,  D2,  Ml,  M2)  all  show  ct  ~  1  at  0.5  K,  slowly  decreasing 
with  temperature  T  to  about  a  ~  0.9  at  20  mK,  while  above  0.7 
K  a  falls  rapidly  as  T  increases. 

In  contrast,  the  smaller  devices  show  a  markedly  lower  than 
1  below  1  K.  Despite  large  differences  in  shape  and  dimensions, 
devices  Jl,  FI,  Kl  and  PI  behave  quite  similarly:  a  ~  0.8 
at  0.5  K,  decreasing  to  0.7  below  100  mK.  Device  C2  shows 
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Fig.  6.  Effective  magnetic  loss  tangent  tan (fim)  at  1  Hz  versus  temperature  T 
for  data  shown  in  Fig.  5.  The  symbols  are  for  devices  D2  -Br-Ml  +-D1  A. 
M2  — D,-KT  s-.-Fl----£,-Gl  ^rand  C2 


T(K) 


Fig.  7.  Slope  a  of  the  l/f“  flux  noise  versus  temperature  for  10  devices.  The 
symbols  are  for  devices  D2  -*-Ml  -ArM2  -*r-Pl  --■E},-K1  — 0-yJL 

---e-,-Fl-  X ,  and  C2  * . 

a  singular  behavior,  with  considerable  variation  around  0.1  K. 
The  device  with  the  smallest  center  hole  (Ol)  showed  a  consis¬ 
tently  low  slope,  reaching  a  minimum  of  about  0.45  at  0.4  K, 
and  increasing  to  ~0.6  at  30  mK.  We  note  that  while  FI  has 
the  smallest  loop  linewidth  (2  /im)  and  an  intermediate  sized 
center  hole  (390  /nil)  it  behaves  more  like  devices  Jl,  K1  and 
PI  than  the  device  with  the  smallest  center  hole  (Ol).  This  sug¬ 
gests  that  the  inner  length  of  the  loop  may  be  an  important  factor 
in  determining  the  slope  of  the  noise.  As  one  sees  from  Table  I, 
however,  an  overriding  distinction  may  be  that  Dl,  D2,  Ml,  M2 
with  a  >  0.92  all  have  large  linewidths  (w  >  300  //in)  while 
the  remaining  devices  with  a  <  0.74  have  relatively  narrow 
linewidths  (w  >  50  /um). 

VII.  Conclusion 

In  conclusion,  we  measured  1  / fv>  flux  noise  in  ten 
Nb  —  NbOx  —  Pbln  dc  SQUIDs  below  1  K  and  found  a  strong 
dependence  of  the  slope  a  on  temperature  and  SQUID  geom¬ 
etry.  Smaller  SQUIDs — with  narrower  linewidths — showed 
markedly  lower  values  of  a  and  a  different  dependence  on  tem¬ 
perature  than  larger  SQUIDs — with  wider  linewidths.  These 
results  are  inconsistent  with  a  random  distribution  of  uncorre- 
lated  noise  sources  [7],  which  would  cause  the  noise  magnitude 
but  not  the  slope  to  vary  with  geometry,  and  suggests  that  there 
is  a  spatial  correlation  to  the  flux  noise.  Indeed,  correlations 
among  spins  are  predicted  to  produce  values  of  a  that  can 
be  well  below  unity  [11],  We  note  that  none  of  our  devices 
showed  the  flattening  of  the  1/f  noise  at  frequencies  below 


fw  predicted  by  the  spin-diffusion  model  [9],  These  authors 
suggest,  however,  that  weakly  coupled  spins  further  away  from 
the  superconductor  and  the  presence  of  strongly  coupled  spin 
pairs  could  produce  noise  with  nonzero  slope  at  frequencies 
below  fw. 

Our  measurements  were  carried  out  long  before  the  develop¬ 
ment  of  microscopic  models  of  the  1/f  flux  noise  [4],  [6],  [9], 
[10],  [11],  Lacking  insight  from  these  models,  we  did  not  sys¬ 
tematically  vary  only  d  and  w,  the  key  parameters  in  the  models 
[9],  [12],  as  was  done,  for  example,  in  [13],  Thus,  we  cannot 
usefully  compare  the  noise  magnitude  of  our  own  devices  with 
predictions  [4],  [9],  [12],  Future  measurements  should  be  car¬ 
ried  out  with  SQUIDs  in  which  only  a  single  parameter  is  varied 
to  determine  whether  a  depends  on,  for  example,  linewidth  or 
hole  size. 
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